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We measured Newton's gravitational constant G using a new torsion balance method. Our 
technique greatly reduces several sources of uncertainty compared to previous measurements: (1) it 
is insensitive to anelastic torsion fiber properties; (2) a flat plate pendulum minimizes the sensitivity 
due to the pendulum density distribution; (3) continuous attractor rotation reduces background 



noise. We obtain G = (6.674215 ± 0.000092) x 10~ n m 3 kg -1 s~' A ; the Earth's mass is, therefore, 
(5.972245 ±0.000082) x 10 24 kg and the Sun's mass is M = (1.988435 ±0.000027) x 10 30 kg. 
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PACS numbers: 06.20.Jr, 04.80.-y, 06.30.Gv, 97.10.Nf, 91.10 



The gravitational constant G along with Planck's con- 
stant h and the speed of light c, is one of the most fun- 
damental and universal constants in nature. Unlike most 
other physical constants, the value of G is not precisely 
known due to the weakness and non-shieldability of grav- 
ity. Since the first laboratory measurement by Cavendish 
over 200 years ago, the reduction in uncertainty in G has 
been only a factor of about ten per century. To make 
matters worse, measurements over the last ten years 
all have larger error-bars than the value of Luther and 
Towler || on which the 1986 accepted value [Q was based. 
In addition, some of these recent measurements lie far 
outside the ±128 ppm uncertainty of the 1986 accepted 
value, bringing this accuracy of G into question. Recog- 
nizing this situation, the currently recommended value is 



G = (6.673±0.010) x 10" n m 3 kg" 



correspond- 



ing to an uncertainty of ±1500 ppm. 

We developed a new technique ^ that offers several 
substantial advantages over previous methods. At the 
heart of the apparatus is a torsion balance placed on 
a turntable located between a set of attractor spheres. 
The turntable is first rotated at a constant rate so that 
the pendulum experiences a sinusoidal torque due to the 
gravitational interaction with the attractor masses. A 
feedback is then turned on that changes the rotation rate 
so as to minimize the torsion fiber twist. The resulting 
angular acceleration of the turntable, which is now equal 
to the gravitational angular acceleration of the pendu- 
lum |jj , is determined from the second time-derivative of 
the turntable angle readout. Since the torsion fiber does 
not experience any appreciable deflection, this technique 
is independent of many torsion fiber properties includ- 
ing anelasticity, which may have led to a bias |||||]] hi 
previous measurements. 

The gravitational angular acceleration a of a torsion 
pendulum expanded in spherical multipole moments is: 
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and the Q/ m multipole fields |T(| of the external mass 
distribution, (j) is the turntable angle with respect to 
the mass distribution, and I is the pendulum moment 
of inertia about the torsion fiber. Since the series con- 
verges rapidly, the biggest term in Eq. [l] is the one with 
I = m = 2: 
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For an ideal, infinitely-thin two-dimensional vertical 
plate the quotient (722 /I is independent of the pendulum 
mass distribution, p(r p ), i.e. mass, size, shape etc.: 
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The flat pendulum geometry presents a significant ad- 
vance over previous torsion balance measurements, where 
the biggest contribution to the systematic uncertainty 
was due to the pendulum's mass distribution. A small 
and easily calculable modification can be applied to the 
ideal ratio for a thin rectangular plate with finite thick- 
ness t and width w: 



q22 = w z -t 1 15 
I ~ w 2 +t 2 \l 32tt 
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where qi m are the multipole moments of the pendulum 



All I > 2, to = 2 terms in Eq. [l] also contribute 
to the signal sin(20). However, these contributions 
are small since the expansion in Eq. [l] converges as 
(R P end/Rattr) l ~ 2 , where R pen d and R a ttr are the char- 
acteristic radii of the pendulum and the source. In addi- 
tion, all I = odd pendulum moments and attractor fields 
are removed by designing the pendulum "up-down" sym- 
metric. The I = 4, to = 2 contribution is also eliminated 
by design: the 942 of a rectangular plate pendulum with 
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height h vanishes by choosing h 2 = jq(w 2 + t 2 ), and 
the Qi2 vanishes by using two spheres on each side with 
their centers spaced vertically by z — ^2/3p, where p 
is the radial distance from the pendulum axis. The low- 
est I > 2, m = 2 contribution which does not vanish by 
design is the small and easily calculable cte2- 

a 62 _ 99 213(w 4 + t 4 ) + 626w 2 t 2 

~ 7683200 / ' ^' 

Table I contains the numeric values of the «62 and as2 
corrections. We checked the accuracy of our multipole 
analysis with a full angular acceleration calculation us- 
ing direct numeric integration. 

The attractor spheres are located on a separate coaxial 
turntable which is rotated with angular velocity to a (t) = 
lu c [ + u>i(t), where u>i{t) is the angular velocity of the 
torsion balance turntable. The difference angular veloc- 
ity tud — <f>, is held constant. Rotation of the attractor 
masses allows us to cleanly remove gravitational interac- 
tions due to the environment. Furthermore, we are able 
to set the signal, sin(2wdi), at a relatively high frequency 
to suppress the 1 //-noise characteristic of the torsion bal- 
ance and the gravitational background. 

A schematic of the apparatus is shown in Fig. 1. The 
torsion balance turntable consisted of an air-bearing, a 
precision angle encoder, and an eddy-current motor. The 
torsion pendulum was located in an aluminum vacuum 
chamber and was surrounded by a /i-metal shield. The 
pendulum was hung from a 41.5 cm long, 17 pm diameter 
tungsten-fiber, which was attached to a swing damper. 
The pendulum was a 1.506 mm thick, 76 mm wide, and 
41.6 mm tall Pyrex-glass plate with a thin gold-coating. 
The small pendulum deflection angle was sensed with an 
autocollimator using four reflections off the pendulum 
plate. 

The centers of the spheres were located at p=16.76 cm 
on three stainless steel seats that were embedded in two 
cast-aluminum shelves. The shelves were supported by 
a turntable made from a steel bearing. The attractor 
spheres were machined from the same selected stock of 
ultrasonically tested #316 stainless steel. Their average 
diameter was 124.89 mm and their mass was ?»8.140 kg 
jnj . A pressure-dependent air-density correction was ap- 
plied to the G-measurements. The apparatus was located 
in the former cyclotron cave in the Nuclear Physics Lab 
at the University of Washington on a massive platform 
3.5 m above the floor. The partially underground room 
was temperature stabilized. Temperature drifts and fluc- 
tuations were typically < 0.05 K/day. The instrument 
itself was in a passive thermal enclosure made of Sty- 
rofoam. The apparatus-temperature was monitored and 
a correction was applied to compensate for the thermal 
expansion of the attractor mass assembly during the mea- 
surement. 

A digital signal processor (DSP) recorded the data and 
controlled the experiment. The DSP's loop frequency 



(2.5 KHz) was used for the timing, which was derived 
from a temperature-controlled quartz oscillator which 
was calibrated with a GPS receiver. The data were aver- 
aged by the DSP over exactly one second and uploaded 
to the host PC. 

We recorded six data runs, each approximately three 
days long. A typical data segment is shown in Fig 2. 
After each data run the spheres were moved to loca- 
tions differing by 90° in azimuthal angle on the attrac- 
tor turntable shelves. In addition, the spheres from the 
upper shelf were placed on the lower shelf and the ori- 
entation of the spheres were changed in order to average 
over density fluctuations and non-sphericity. The posi- 
tion of the spheres were measured before and after each 
data run. We then repeated the entire measurement cy- 
cle using four different spheres. 

Our largest systematic uncertainty was due to the at- 
tractor mass distance measurement. We used a spe- 
cially fabricated micrometer tool made primarily of In- 
var for the horizontal distance measurement between the 
spheres. Before and after each distance measurement the 
tool was compared to an Invar ball bar standard that was 
calibrated at NIST to within 0.2 ^m Q. The vertical 
spacing between the spheres was measured by inserting a 
small gauge block that was ssl0-20 /im thinner than the 
gap between the sphere surfaces. We inferred the spacing 
by optically measuring the angle through which the gauge 
block could be tilted. The temperature of the attractor 
mass was recorded during the distance measurements. 

For most of the data u>i was set to «5.3 mrad/s and Ud 
was set to 20.01015 mrad/s so that the signal frequency 
occurred at ~6.37 mHz. The rotation frequencies were 
chosen to be incommensurate. We tested a wide variety 
of LOi and u>d and found the results independent of these 
angular velocities. 

The lab-fixed horizontal magnetic field at the appa- 
ratus was measured to be «100 mG. We ran tests by 
exaggerating the field at the center of the apparatus to 
«5 G and at the location of the spheres to «100 G. The 
observed acceleration difference due to the exaggerated 
fields was (6 ± 8) x 10~ 12 rad/s 2 . Therefore, a 0.6 ppm 
error was attributed to magnetic accelerations. 

To investigate the rotating temperature gradient sen- 
sitivity, two 5W-heaters were installed on the turntable. 
These rotating heaters exaggerated the normally ob- 
served temperature variation at 2u>d by greater than a 
factor of 250. The acceleration difference amplitude with 
the heaters activated was | Aa| = (21±22) x 10 -12 rad/s 2 . 
Therefore the rotating temperature gradient coupled ac- 
celeration was < 2 x 10 -13 rad/s 2 (0.4 ppm). 

The torsion balance turntable angle was numerically 
differentiated twice to yield angular acceleration. All the 
data were divided into segments comprised of twenty si- 
nusoidal cycles and fitted using a least squares method. 
The fitting function included the signal and its harmon- 
ics, the room background and its harmonics, an offset 
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and a linear drift. The statistical error was derived from 
the scatter of the individual fit values. 

Tabic 1 lists the corrections that were applied to the 
data. Since the DSP averaged over r = 1 s intervals 
the true amplitudes are attenuated by a small amount: 
Sm i"r T ^ • Also the numerical derivative of sm(2ujdt) is 
smaller than the analytic derivative by Sln ^"^^ , where 
At = 10 s was the chosen time increment. Since this 
correction was substantial we tested a wide variety of in- 
crements and found that our results were independent of 
the choice of At. We verified the accuracy of the data 
analysis by simulating the data numerically using var- 
ious levels of drifts, l//-noise, and gravitational room 
background. 

The overall linearity and insensitivity to the signals at 
Auid was tested using two pairs of spheres separated by 
45° on each side. These 8-sphere results agreed with the 
regular 4-sphere results within the statistical errors. 

Due to the finite gain of the acceleration feedback loop 
a small twist remained in the fiber. The twist angle at 
2u>d was essentially unresolved so that a 0.35 ppm error 
was used instead of a correction. 

We combined the results from the different attractor 
mass configurations into pairs that averaged out acceler- 
ations due to the attractor turntable itself and the com- 
bination of three pairs yielded our G- value (Fig. 3). Af- 
ter combining the two G-values obtained with different 
spheres our value for the gravitational constant is 



G = (6.674215 ± 0.000092) x 10 



-li 



m 3 kg 1 s 2 . 



(7) 



The systematic and statistical uncertainties are shown in 
Table 2. The total uncertainty is the quadrature sum of 
the individual uncertainties. Fig. 4 compares our results 
to other recent measurements. 

Combining our value of G with the geocentric value 
GM§ we determine the mass of the Earth to be 

M e = (5.972245 ± 0.000082) x 10 24 kg. (8) 

Likewise, using GMq the mass of the Sun is 

Mq = (1.988435 ± 0.00 2 7) x 10 30 kg. (9) 
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TABLE I. Summary of correction factors. 



finite pendulum thickness 1.0007857 

pendulum attachment and imperfections 1.0000433 

^-correction 0.9998767 

a 8 2-correction 0.9999951 

data averaging (r = 1 s, u>d = 20 mrad/s) 1.0000667 
numeric derivatives (At = 10 s, cud = 20 mrad/s) 1.0134544 

Total: 1.0142322 
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TABLE II. One a error budget. 



quantity 


measurement 
uncertainty 


AG/G 


(ppm) 


Systematic errors: 
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A n / im 
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attractor masses: 
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ball-bar calibration 


^ n o i,™ 
<. U.^i /im 


1 A 
1.4 


vertical separation 


< 1.0 /im 


5.2 


sphere diameter 


< 1.5 /im 


2.6 


temperature uncertainty 


< 100 mK 


6.9 


mass 


< 3.0 mg 


0.4 


air humidity 




0.5 


residual twist angle: 




0.3 


magnetic fields: 




0.6 


rot. temperature grad.: 




0.4 


time base: 


< 10" 7 


0.1 


data reduction: 




2.0 


Statistical error: 




5.8 


Total: 




13.7 
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FIG. 2. Raw acceleration data: a half-hour segment of the 
twice numerically differentiated turntable angle. The signal 
frequency was constant and could be freely selected. 
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FIG. 3. The results of two data sets taken with different 
spheres. Each data point is the combination of a pair of at- 
tractor configurations that together eliminate accelerations 
due to the attractor turntable. Combining the three pairs 
taken with different sphere orientations optimally reduced ef- 
fects from sphere-density and shape imperfections. The dis- 
played uncertainties are statistical only. 
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FIG. 4. Comparison to other measurements 
after 1995 and with AG/G < 1000 ppm. 
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FIG. 1. Cut-away view of the apparatus. 
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